Distinct organization of adaptive immunity in the long-lived rodent Spalax galili.

A balanced immune response is a cornerstone of healthy aging. Here, we uncover distinctive features of the long-lived blind mole-rat (Spalax spp.) adaptive immune system, relative to humans and mice. The T-cell repertoire remains diverse throughout the Spalax lifespan, suggesting a paucity of large long-lived clones of effector-memory T cells. Expression of master transcription factors of T-cell differentiation, as well as checkpoint and cytotoxicity genes, remains low as Spalax ages. The thymus shrinks as in mice and humans, while interleukin-7 and interleukin-7 receptor expression remains high, potentially reflecting the sustained homeostasis of naive T cells. With aging, immunoglobulin hypermutation level does not increase and the immunoglobulin-M repertoire remains diverse, suggesting shorter B-cell memory and sustained homeostasis of innate-like B cells. The Spalax adaptive immune system thus appears biased towards sustained functional and receptor diversity over specialized, long-lived effector-memory clones-a unique organizational strategy that potentially underlies this animal's extraordinary longevity and healthy aging.

[A Pipeline for the Error-free Identification of Somatic Alu Insertions in High-throughput Sequencing Data].

Retroelements are considered as one of the important sources of genomic variability in modern humans. It is known that transposition activity of retroelements in germline cells generates new insertions in various genomic loci and sometimes results in genetic diseases. Retroelements activity in somatic cells is restricted by different cellular mechanisms; however, there is an evidence for it in some tissue types. Somatic insertions can trigger tumorigenesis or participate in normal functioning such as generation of neurons' plasticity. In spite of the rapid development of high-throughput sequencing methods a confident detection of somatic insertions is still quite a challenging task. That, in part, is due to the absence of adequate bioinformatic tools for the analysis of sequencing data. Here, we propose an advanced computational pipeline for the identification of somatic insertions in datasets generated by selective amplification and high-throughput sequencing of genomic regions flanking insertions of AluYa5. Particular attention is paid for the identification of various artifacts arising in course of library preparation and the parameters for their filtration. Pipeline sensitivity is confirmed by in silico experiments with artificial datasets. Using the proposed pipeline we remove at least 80% of artifacts and preserve 75% of potentially somatic insertions. The approaches used in this work can be applied for the study of other mobile elements insertion variability.

Characterization of the T-cell Repertoire after Autologous HSCT in Patients with Ankylosing Spondylitis.

Autologous hematopoietic stem cell transplantation (HSCT), a safer type of HSCT than allogeneic HSCT, is a promising therapy for patients with severe autoimmune diseases (ADs). Despite the long history of medical practice, structural changes in the adaptive immune system as a result of autologous HSCT in patients with various types of ADs remain poorly understood. In this study, we used high-throughput sequencing to investigate the structural changes in the peripheral blood T-cell repertoire in adult patients with ankylosing spondylitis (AS) during two years after autologous HSCT. The implementation of unique molecular identifiers allowed us to substantially reduce the impact of the biases occurring during the preparation of libraries, to carry out a comparative analysis of the various properties of the T-cell repertoire between different time points, and to track the dynamics of both distinct T-cell clonotypes and T-cell subpopulations. In the first year of the reconstitution, clonal diversity of the T-cell repertoire remained lower than the initial one in both patients. During the second year after HSCT, clonal diversity continued to increase and reached a normal value in one of the patients. The increase in the diversity was associated with the emergence of a large number of low-frequency clonotypes, which were not identified before HSCT. Efficiency of clonotypes detection after HSCT was dependent on their abundance in the initial repertoire. Almost all of the 100 most abundant clonotypes observed before HSCT were detected 2 years after transplantation and remained highly abundant irrespective of their CD4+ or CD8+ phenotype. A total of up to 25% of peripheral blood T cells 2 years after HSCT were represented by clonotypes from the initial repertoire.

Tracking T-cell immune reconstitution after TCRαß/CD19-depleted hematopoietic cells transplantation in children.

αßT-cell-depleted allogeneic hematopoietic cell transplantation holds promise for the safe and accessible therapy of both malignant and non-malignant blood disorders. Here we employed molecular barcoding normalized T-cell receptor (TCR) profiling to quantitatively track T-cell immune reconstitution after TCRαß-/CD19-depleted transplantation in children. We demonstrate that seemingly early reconstitution of αßT-cell counts 2 months after transplantation is based on only several hundred rapidly expanded clones originating from non-depleted graft cells. In further months, frequency of these hyperexpanded clones declines, and after 1 year the observed T-cell counts and TCRß diversity are mostly provided by the newly produced T cells. We also demonstrate that high TCRß diversity at day 60 observed for some of the patients is determined by recipient T cells and intrathymic progenitors that survived conditioning regimen. Our results indicate that further efforts on optimization of TCRαß-/CD19-depleted transplantation protocols should be directed toward providing more efficient T-cell defense in the first months after transplantation.

High-quality full-length immunoglobulin profiling with unique molecular barcoding.

High-throughput sequencing analysis of hypermutating immunoglobulin (IG) repertoires remains a challenging task. Here we present a robust protocol for the full-length profiling of human and mouse IG repertoires. This protocol uses unique molecular identifiers (UMIs) introduced in the course of cDNA synthesis to control bottlenecks and to eliminate PCR and sequencing errors. Using asymmetric 400+100-nt paired-end Illumina sequencing and UMI-based assembly with the new version of the MIGEC software, the protocol allows up to 750-nt lengths to be sequenced in an almost error-free manner. This sequencing approach should also be applicable to various tasks beyond immune repertoire studies. In IG profiling, the achieved length of high-quality sequence covers the variable region of even the longest chains, along with the fragment of a constant region carrying information on the antibody isotype. The whole protocol, including preparation of cells and libraries, sequencing and data analysis, takes 5 to 6 d.

Advanced lymphoblastic clones detection in T-cell leukemia.

T cell acute lymphoblastic leukemia (T-ALL) is an aggressive malignant neoplasm of the lymphocyte precursors that suffered malignant transformation arresting the lymphoid cell differentiation. Clinical studies revealed monoor, more rarely, oligoclonal nature of the disease. A precise identification of malignant clone markers is both the crucial stage of early diagnostics and the essential prognostic factor for therapeutic treatment. Here we present an improved system for unbiased detection of lymphoblastic clones in bone marrow aspirates of T-ALL patients. The system based on multiplex PCR of rearranged T-cell receptor locus (TRB) and straightforward sequencing of the resulted PCR fragments. Testing of the system on genomic DNA from Jurkat cell line and four clinical bone marrow aspirates revealed a set of unique TRB rearrangements that precisely characterize each of tested samples. Therefore, the outcome of the system produces highly informative molecular genetic markers for further monitoring of minimal residual disease in T-ALL patients.

[A novel approach to identification of somatic retroelements' insertions in human genome].

The activity of retroelements is one of the factors leading to genetic variability of the modern humans. Insertions of retroelements may result in alteration of gene expression and functional diversity between cells. In recent years an increasing amount of data indicating an elevated level of retroelements' mobilisation in some human and animal tissues has been reported. Therefore, the development of a system for the detection of somatic retroposition events is required. Here we describe a novel approach to the whole-genome identification of somatic retroelemts' insertions in human genome. The developed approach was applied for the comparisons of somatic mosaicism levels in two tissues of the investigated individual. A total of 3410 insertions of retroelements belonging to AluYa5 subfamily were identified.

[The oligoclonal expansion of T cells: study of its stability over time].

A novel experimental approach to the investigation of the repertoire of peripheral T lymphocytes of patients suffering from ankylosing spondylitis (AS) is proposed. This approach is based on the wide-range sequencing of cDNA of the beta-chain of the T-cellular receptor (TcR). The results of the analysis of the diversity of sequences of the TcR antigen-binding domain (CDR3) inside the total pool of one patient with AS are presented by the example of the second V family (BV2) of TcR. The expansion of six independent TcR-expressing clones of T cells with a similar amino acid sequence of the CDR3 domains was proposed based on the results of the comparative structural analysis of the clone libraries of the cDNA of TcR BV2. The long-time stable expansion of these T clones was demonstrated during the development of the disease by specific monitoring.

Association of ERAP1 Allelic Variants with Risk of Ankylosing Spondylitis.

Ankylosing spondylitis (AS) belongs to a group of autoimmune diseases affecting the axial skeleton. Beside thehla-b*27allele, several other human genes that control the variety processes of immune homeostasis are considered to be associated with AS manifestation in different human populations. Among strong associated non-MHC geneserap1 encodingthe endoplasmic reticulum aminopeptidase 1 isoform was recently identified by single nucleotide polymorphisms (SNPs) meta analysis. In our study we inspected the genetic association of five non-synonymous coding SNPs fromerap1 withAS in Caucasians. We implemented the SSP-PCR system for precise genotyping of 87hla-b*27positive AS patients and 77hla-b*27healthy donors from the Russian population. Considerable differences in allele's frequencies within patients vs control cohort were shown for 3 of 5 SNPs under investigation. Using the EM-algorhitm we reconstructed 3-marker haplotypes that distinguish with high probability two cohorts due to differences in the haplotypes frequencies. In such a way both the sensitive, CCT, haplotype and the protective, TTC, one were predicted. To verify the calculation we determined genuine frequencies of 5-marker haplotypes in AS cohort by haplotyping of individual cDNA samples using improved SSP-PCR primer set. We demonstrated that the frequencies ofin silicareconstucted haplotypes and the frequencies of experimentally detected haplotypes are in a good agreement. Frequency of the risk haplotype CCT (rs17482078/10050860/2287987) detected within AS cohort reaches 88%, as well as the frequency calculated by EM-algorhitm.

Individual characterization of stably expanded T cell clones in ankylosing spondylitis patients.

Ankylosing spondylitis (AS) is commonly characterized by clonal expansions of T cells. However, these clonal populations are poorly studied and their role in disease initiation and progression remains unclear. Here, we performed mass sequencing of TCR V beta libraries to search for the expanded T cell clones for two AS patients. A number of clones comprising more than 5% of the corresponding TCR V beta family were identified in both patients. For the first time, expanded clones were shown to be stably abundant in blood samples of AS patients for the prolonged period (1.5 and 2.5 years for two patients, correspondingly). These clones were individually characterized in respect to their differentiation status using fluorescent cell sorting with CD27, CD28, and CD45RA markers followed by quantitative identification of each clone within corresponding fraction using real time PCR analysis. Stable clones differed in phenotype and several were shown to belong to the proinflammatory CD27 - /CD28 - population. Their potentially cytotoxic status was confirmed by staining with perforin-specific antibodies. Search for the TCR V beta CRD3 sequences homologous to the identified clones revealed close matches with the previously reported T cell clones from AS and reactive arthritis patients, thus supporting their role in the disease and proposing consensus TCR V beta CDR3 motifs for AS. Interestingly, these motifs were also found to have homology with earlier reported virus-specific CDR3 variants, indicating that viral infections could play role in development of AS.

[Retroposons in modern human genome evolution].

The ascertainment of the rates and driving forces of human genome evolution along with the genetic diversity of populations or separate population groups remains a topical problem of fundamental and applied genomics. According to the results of comparative analysis, the most numerous human genome structure peculiarities are connected with the distribution of mobile genetic retroelements - LTR, LINE1, SVA, and Alu repeats. Due to the wide distribution in different genome loci, conversed retropositional activity, and the retroelements regulatory potential, let us regard them as one of the significant evolutionary driving forces and the source of human genome variability. In the current review, we summarize published data and recent results of our research aimed at the analysis of the evolutionary impact of the young retroelements group on the function and variability of the human genome. We examine modern approaches of the polygenomic identification of polymorphic retroelements inserts. Using an original Internet resource, we analyze special features of the genomic polymorphic inserts of Alu repeats. We thoroughly characterize the strategy of large-scale functional analysis of polymorphic retroelement inserts. The presented results confirm the hypothesis of the roles of retroelements as active cis regulatory elements that are able to modulate surrounding genes.

[A new database on polymorphic retroelements in human genome (PRED)].

Comparison of primate genomes sequences has confirmed the evidence that substantial part of intra- and interspecies differences is provided by retroelements. Human genome contains thousands of polymorphic retroelement copies considered to be perspective molecular genetic markers of new generation. However utilization of polymorphic retroelements as molecular genetic markers is limited due to lack of systematic data on their number, genomic context and distribution among human populations. We have created first bilingual (Russian/English) internet-resource devoted to known polymorphic retroelements discovered in human genome by our group as well as by other researchers worldwide. The database contains information about each retroelement copy location, position relative to known and predicted genes, frequency of alleles in human populations and others. Our internet portal allows to perform a search in database using multiple search conditions and available on http://labcfg.ibch.ru/home.html. The database provides an opportunity to investigate distribution of polymorphic retroelements in human genome and to design new genetic markers for various population and medical studies.

Full-sized HERV-K (HML-2) human endogenous retroviral LTR sequences on human chromosome 21: map locations and evolutionary history.

One of the evolutionary mechanisms for acquisition of novel functional sequences can be domestication of exogenous retroviruses that have been integrated into the germ line. The whole genome mapping of such elements in various species could reveal differences in positions of the retroviral integration and suggest possible roles of these differences in speciation. Here, we describe the number, locations and sequence features of the human endogenous retrovirus HERV-K (HML-2) long terminal repeat (LTR) sequences on human chromosome 21. We show that their distribution along the chromosome is not only non-random but also roughly correlated with the gene density. Amplification of orthologous LTR sites from a number of primate genomes produced patterns of presence and absence for each LTR sequence and allowed determination of the phylogenetic ages and evolutionary order of appearance of individual LTRs. The identity level and phylogenetic age of the LTRs did not correlate with their map locations. Thus, despite the non-random distribution of LTRs, they have apparently been inserted randomly into the chromosome relative to each other. As evidenced in previous studies of chromosomes 19 and 22, this is a characteristic of HERV-K integration.